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We have measured temperature dependence of low-frequency Raman spectra in 



CuGe03, and have observed the quasi-elastic scattering in the (c, c) polarization 
above the spin-Peierls transition temperature. We attribute it to the fluctuations of 
energy density in the spin system. The magnetic specific heat and an inverse of the 
magnetic correlation length can be derived from the quasi-elastic scattering. The 
inverse of the magnetic correlation length is proportional to (T — Tsp) 1 ^ 2 at high 
temperatures. We compare the specific heat with a competing- J model. This model 
cannot explain quantitatively both the specific heat and the magnetic susceptibility 
with the same parameters. The origin of this discrepancy is discussed. 

PACS numbers: 78.30.-j, 78.20.-e, 75.90,+w 



1. INTRODUCTION 

Hase, Terasaki, and UchinokuraS reported that the magnetic susceptibility in all directions 
in CuGe03 decreases exponentially below T$p = 14 K, and they concluded that it is due to 
a spin-Peierls (SP) transition. In the SP phase, an energy gap between the singlet ground 
and the triplet excited states opens and a lattice dimerization is formed.i After this report, 
many experimental studies have been performed in this compound to reveal the nature of 
the SP phase. The magnetic excitations and the SP gap have been extensively studied by 
inelastic-neutroni and Raman scattering.^ Moreover the lattice dimerization was observed 
below T,sp,!~i and then the SP transition was established in CuGeOs. 

On the other hand, Hase, Terasaki, and Uchinokura0 pointed out that the magnetic sus- 
ceptibility above T$p cannot be described quantitatively by Bonner-Fisher's theory! based 
on an S = 1/2 one-dimensional (ID) Heisenberg antiferromagnetic (AF) Hamiltonian with 
a nearest-neighbor (nn) exchange interaction. Recently, Riera and Dobryi and Castilla, 
Chakravarty, and Emery0 proposed independently for CuGeO"3 an S = 1/2 ID Heisenberg 
model with the nearest- and next-nearest-neighbor (nnn) AF exchange interactions in the 
chain (a competing- J model). The Hamiltonian of this model is written as 

N 

H = J (Si • S i+ i + aSi ■ S i+2 ) , (1) 

i=l 

where J (> 0) and a (> 0) are the nn AF exchange interaction and the ratio of the nnn 
interaction to the nn one, respectively. In this model, if a is larger than a critical value a c , 
an energy gap between the singlet ground and the triplet excited states appears without any 
lattice dimerization. a c is estimated to be O.2411.0'0 Riera and Dobry suggested that the 
model with J = 160 K and a = 0.36 can explain the magnetic susceptibility above T^p.i On 
the other hand, Castilla, Chakrabarty, and Emery reported that the model with J = 150 K 
and a = 0.24 can explain both the magnetic susceptibility above T$p and a dispersion of 
magnetic excitations at low temperatures.0 

In order to check the validity of this model, other experimental results will be compared 
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with this model. As is well known, the magnetic specific heat (C m ), as well as the mag- 
netic susceptibility, is an important physical quantity Therefore, it is necessary to measure 
accurately C m at high temperatures. However, C m at high temperatures has been poorly 
studied, although the anomalous specific heat around Tsp in CuGeC>3 was investigated in 
many experiments. EUll Weiden et a/.@ observed the specific heat below 300 K in poly- 
crystalline CuGeC>3 using a continuous-heating adiabatic method. They obtained C m after 
subtracting a calculated lattice specific heat from the observed one. The lattice specific heat 
was assumed to be written as a sum of Debye functions. They reported that C m shows 
a broad maximum around 40 K. As we will describe below, our result contradicts theirs. 
As is well known, C m at high temperatures cannot be obtained accurately by this method 
because effects of optical phonons, which play an important role in the specific heat at high 
temperatures, were not taken into account precisely. 

It is noted that the quasi-elastic scattering is a powerful tool to derive precisely C m 
especially at high temperatures. In our brief report ,0 a large tail near an incident laser line 
with a strong temperature dependence was observed at high temperatures in single- crystal 
CuGeO"3, the profile in the Stokes and anti-Stokes scattering was fit to a Lorentzian curve. 
We assigned it to the quasi-elastic scattering due to spin-energy fluctuations. We can obtain 
roughly the magnetic specific heat from the quasi-elastic Raman scattering in CuGeO"3, 
because the integrated intensity of the quasi-elastic scattering is proportional to C m T 2 where 
T is the temperature.0 The quasi-elastic scattering was applied to some magnets such as a 
two-dimensional antiferromagnet FePS3,il and later a ID antiferromagnet KCuF 3 .il 

In our brief reportS we observed Raman spectra in both the Stokes and anti-Stokes 
regions at seven different temperatures. The Raman spectra, however, contained a weak 
component of the direct scattering, judging from the fact that a plasma line of the incident 
laser line appeared at -13.4 cm -1 . The quality of the sample was probably worse than that 
used in the present work. It is, of course, necessary to use high-quality samples in order to 
the effects of the direct scattering when we measure the quasi-elastic scattering. Then we did 
not obtain accurately the integrated intensity of the quasi-elastic scattering, in particular at 
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low temperatures. At low temperatures the intensity of the quasi-elastic scattering decreased 
in CuGe0 3 and it was considerably affected by the direct scattering, when we estimated the 
integrated intensity in our brief report. 

Very recently, on the other hand, van Loosdreht et ai.@ measured Stokes Raman spec- 
tra with a backscattering geometry and proposed that the low-frequency Raman spectrum 
at high temperatures originated from the diffusive-type fluctuations in the four-spin time 
correlation functions. Richards and Brya@ theoretically studied Raman scattering from the 
diffusive-type fluctuations in the spin system by using a four-spin time correlation function 
and obtained a result that the profile was described as a Gaussian-like line at uj=0 at high 
temperatures. Therefore the model proposed by van Loosdreht et at seems to contradict 
the interpretation for the quasi-elastic scattering in our brief report .0 

Therefore the careful and precise measurement of the low-frequency Raman spectrum in 
both the Stokes and anti-Stokes regions is needed in a high-quality sample in order to reveal 
clearly the origin of the quasi-elastic Raman scattering. In the present paper we shall obtain 
the magnetic specific heat precisely above T$p and study effects of the nnn AF exchange 
interaction along the magnetic chain on the magnetic properties of CuGeC>3. 



2. EXPERIMENTAL DETAILS 

A single crystal used in this study was made by a floating-zone method. The dimension 
of the crystal is 1.2 x 1.8 x 10.7 mm 3 parallel to the a, b, and c axes, respectively.!^ The 
single crystal has flat ac and be faces. On the other hand, it is difficult to polish the ab 
face because of cleavage on the be face. However, effects of local strain on the ab face were 
negligible when the light was incident and scattered on the ac and be faces. The be face was 
attached to the sample holder by dilute varnish. 

Raman-scattering experiments were carried out using the 5 145- A line of an Ar + ion 
laser. The polarized quasi-elastic Raman scattering was measured in either a right-angle or 
a quasi-back scattering geometry. The detailed temperature dependence of the quasi-elastic 
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Raman scattering was measured in the right-angle scattering geometry. To eliminate plasma 
lines and to get a complete linear polarization, we made the laser light pass through a filter 
and a polarizer and then the laser light was incident on the sample in a cryostat. The 
light was incident on the etc face for the right-angle scattering and on the be face for the 
quasi-back scattering. The scattered light was dispersed by a Jobin-Yvon U1000 double- 
grating monochromator and was detected by a photon counting system. The whole system 
is controlled by a microcomputer. To measure Raman spectra below room temperature, 
we used a closed-cycle cryostat. The temperature of the sample is controlled by a PID 
temperature controller. The temperature of the sample is equal to that of the sample holder 
within ±0.1 K. Since the sample is transparent, the local temperature rise due to the incident 
light is negligible. 



3. THEORIES 



A. QUASI-ELASTIC SCATTERING 



We summarize briefly the theory of the quasi-elastic light scattering in antiferromagnets. 
According to Halley's theory,!! the intensity of the quasi-elastic light scattering is given by 
the following Fourier component of the correlation function of the magnetic energy density: 

I{u>) = 7 / dte-^(E(q, t)E*(-q, 0)> , (2) 



where E(q,t) and Tiuj denote the magnetic energy density with the scattering vector q 
and the energy transfer, respectively. 7 is the q- and ^-independent but polarization- 
dependent coefficient. We assume that 7 is temperature-independent. (■ ■ •) represents a 
thermal average. In the hydrodynamic condition,!! Eq. (fj) is transformed to 

77 1 = 1UJ C m TD T q 2 

[UJ) 1 - e-^I^T ^2 + (D T q 2 ) 2 • 1 } 

Here is a thermal diffusion constant which is given by 

D T = K/C m , (4) 



where K is the magnetic contribution to the thermal conductivity, q is approximately equal 
to |q | sin 9/2, where qo and 9 are the wave vector of the incident light and the scattering 
angle, respectively. In a high-temperature approximation, Eq. (|3]) can be written as the 
following Lorentzian curve 

Tt \ i C m T 2 D T q 2 

= 7 + (D T q^ ' (5) 

where 7' = ks^/h. Using Eq. (g), we can estimate the magnetic specific heat and the 
thermal diffusion constant from the integrated intensity and the half width at half maxima 
of the quasi-elastic scattering, respectively. 

B. MAGNETIC SPECIFIC HEAT OF SPIN CHAIN WITH COMPETING EX- 
CHANGE INTERACTIONS 

We can calculate the magnetic specific heat using all energy levels in the spin system. Bonner 
and Fisher calculated magnetic susceptibility and magnetic specific heat in the spin system 
whose Hamiltonian is expressed by Eq. ([!]) with a = and 3 < N < ll.i They calculated all 
energy levels under the periodic condition by means of an exact diagonalization. Because a 
finite-size effect appears at low temperatures, physical quantities at iV — > 00 were estimated 
from an extrapolation. The magnetic specific heat of an infinite ID AF chain with the nn 
AF exchange interaction has a broad maximum at fc^T w 0.481J. 

In the case of a = 0.5 (Majumder-Ghosh model), the ground state of the system with 
4 < N < 8 was studied analytically, and was revealed to be products of dimer states with 
double degeneracy.il This model was developed to the case of iV — » 00 and higher sM An 
energy gap between the singlet ground and the triplet excited states opens in this case, while 
the magnetic excitation is gapless for a = 0. Therefore, a phase transition takes place at 
< a < 1/2. Tonegawa and Harada calculated the energy gap between the singlet ground 
and the triplet excited states numerically when iV < 20,0 and after this work, Okamoto and 
Nomura found that a r at iV — > 00 is 0.2411 from their calculation when iV < 20.0 
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We calculated all energy levels by means of the exact diagonalization for various a 
when 4 < iV < 14. We used a standard subroutine of the Householder method for a 
numerical diagonalization. We calculated the magnetic specific heat {C m ^) for N spins and 
investigated Dn = [NC my N + (N — l)C mi jy-i]/(2iV — 1) vs. N. almost converges a 
value (i.e., C m ) up to N = 14 at least when T > 0.5T max . Here T max is the temperature at 
which C m has a maximum value C m ^ max . Therefore we show for N = 14 as C m in Fig. 
la. T max and C m ^ max decrease and the peak becomes broader with increasing a. We also 
calculate C m /C m , max vs. T/T max to obtain J and a from our observed data, and it is shown 
in Fig. lb. 

4. RESULTS AND DISCUSSION 

Figure 2a shows polarized Raman spectra of CuGeO"3 at room temperature between -100 
and 100 cm -1 strongly in the right-angle scattering geometries. The Raman spectra in the 
regions of u > and uj < correspond to the Stokes and the anti-Stokes components, 
respectively. Since the intensity of the direct scattering is very strong, we cannot obtain 
Raman spectra in the region of |u;| < 7 cm -1 . The quasi-elastic scattering is strongly 
observed in the (c, c) polarization, while it is not observed in the offdiagonal geometries, 
i.e., the (a, b), (a,c), and (c, b) ones. In the SP phase, we observed two-magnetic-exciton 
Raman spectrum below 250 cm' 1 strongly in the (c, c) polarization.^ Thus, the magnetic 
fluctuations should appear strongly in the (c, c) polarization above Tsp, and we conclude 
that the quasi-elastic scattering of CuGeC>3 originates from the spin fluctuations in the spin 
system. 

Figure 2b shows polarized Raman spectra between -100 and 500 cm -1 at room tem- 
perature in the right-angle (fe(c, c)a) and the quasi-back (a(c, c)a and a(c, c)a) scattering 
geometries. Because the Raman spectra in the b(c,c)a and the a(c,c)a scattering geome- 
tries were obtained with different scattering volumes, we cannot compare these scattering 
efficiencies directly. Then we normalized these Raman intensities by the integrated intensity 
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of the 181-cm" 1 A g phonon mode.u On the other hand, since the a(b, b)a Raman spectrum 
was measured in the same experimental condition as the a(c, c)a spectrum, we can compare 
these intensities directly. The 181-cm -1 A g phonon mode may have different intensities in 
the a(b, b)a and the a(c, c)a geometries. Two peaks at 114 and 221 cm -1 are B\ g modes, 
which are basically forbidden in this geometry. One can see a larger tail near the laser line 
and a stronger background at high-frequency region in the a(c, c)a geometry when compared 
with those in the b(c, c)a geometry. We also observed a plasma line of the incident laser at 
-13.4 cm -1 in the quasi-back scattering geometry, while we did not observe it in the right- 
angle scattering geometry. These facts obviously indicate that the quasi-elastic scattering 
was strongly influenced by the direct scattering in the quasi-back scattering geometry. In 
the b(c, c)a geometry, one can notice a crossover between the quasi-elastic scattering spectra 
and the incident laser line at ± 7 cm -1 . 

In our previous report,@ we observed the plasma line at -13.4 cm -1 even in the right- 
angle scattering geometry, because the sample probably contains local strains and defects, 
and its quality was worse than that used in the present work. Moreover, since the ab face 
of the sample used in the previous work0 was not as-grown but cut and polished, the 
effects of local strain on this surface and of strong diffused reflections of the laser light 
from this surface could not be avoided completely. Therefore, the Raman spectra in the 
previous report probably contained the weak component of the direct scattering, and we 
could not obtain the accurate integrated intensity of the quasi-elastic scattering especially 
at low temperatures. 

We measured the temperature dependence of the low-frequency Raman spectra in the 
right-angle scattering geometry using the high-quality sample with as-grown ab surfaces. We 
think that we have done all we could do to avoid the effects of the direct scattering in the 
present measurements. 

Figure 3 shows typical spectra of the quasi-elastic scattering at 250, 120, 50, and 20 
K. The intensity of the quasi-elastic scattering drastically decreases with decreasing tem- 
perature, and we did not observe the quasi-elastic scattering below Tsp. Since we did not 
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observe any asymmetric spectra above 20 K, the observed spectra are fitted to the following 
Lorentzian-type spectral function, instead of Eq. (||), using the method of least squares: 
k 2 T 

I{uj) = 9 ■ + background, (6) 
Ctr + 1 z 

where V and k in Eq. are the damping constant and the coupling coefficient, respectively. 
Comparing Eqs. (|5]) and (|6|), C m and Dt can be evaluated from k 2 and T. The background 
is assumed to be expressed as Auo + B. In the fitting, we use the Raman spectra of \oj\ > 9 
cm -1 , and then we can avoid the effects of the incident laser line near u ~ 0. Below 40 K, 
we could not fit the observed data because their intensities were very weak. 

Van Loosdreht et al. proposed that the origin of the low-frequency Raman spectra is 
the diffusive-type fluctuations in the four-spin time correlation functions.0 According to 
Richards and Bryail, this leads to a Gaussian-like low-frequency Raman spectrum. Then, 
we tried to fit the observed spectra to the following Gaussian-type spectral function: 



k% ( u 2 

2 



I ^ = /F^-n exp I _ ^2" ) + background, (7) 



where kc and Tq corresponds to k and T in Eq. (||). As seen clearly in Fig. 3, the 
Lorentzian-type spectral functions could reproduce the observed data much better than the 
Gaussian-type ones. Especially, the observed spectra near the incident laser line in the 
frequency region of < 15 cm -1 did not fit to the Gaussian-type spectral function. This 
result indicates that the observed quasi-elastic Raman scattering does not originate from 
the diffusive-type four-spin fluctuations but from the spin-energy fluctuations. 

We show the temperature dependence of T and k 2 above 50 K in Figs. 4a and 4b. As is 
drawn by a solid curve in Fig. 4a, T obeys (T-T SP ) 1/2 , and the same result was obtained in 
our brief report. It is noted that V is proportional to Dt, and therefore to an inverse of a 
magnetic correlation length (T oc Dt oc Although we did not observe the quasi- 

elastic scattering below 40 K, T seems to become zero around T$p. Near the critical point, 
the thermal conductivity K which depends on short-range behavior of the system remains 
finite, while the magnetic specific heat C m diverges.il Then D T {= K/C m ) is predicted to 
vanish at T$p. 



Pouget et alu observed the structural fluctuations below about 40 K by means of X-ray 
diffuse scattering and they reported that the inverse correlation length of the structural 
fluctuations was proportional to (T — Tsp) 1 ^ 2 . Recent neutron-scattering studies^S showed 
that the magnetic correlation length agrees with the structural one near T$p. The present 
result indicates that the obtained T of the quasi-elastic scattering above 50 K has the same 
temperature dependence. Then the inverse magnetic correlation length £ _1 is described as 
(T — Tsp) 1 / 2 not only just above Tsp but also at high temperatures. 

Figure 5 shows k 2 /T 2 above 50 K, which is proportional to C m . As is seen in this 
figure, the data have some errors. Thus, we make smoothed data by a five point weighted 
mean (a solid curve), and compare them with theories. These errors are mainly caused 
by the difficulty in the temperature-dependence measurements. In the Raman-scattering 
measurements, it is not so easy to keep the exact alignment, because the spot of the incident 
laser line on the sample is moved by a thermal expansion of the cryostat. 

A broad maximum is seen in the smoothed curve around 90 K. In Fig. 5, we show two 
theoretical magnetic specific-heat curves with J = 237 K and a = 0.30 and with J = 261 
K and a = 0.40, which have maxima at 90 K. These curves are normalized in order that 
the maxima {C m ^ max ) of both the experimental and theoretical C m agree with each other. 
Since the smoothed data exist between two theoretical C m , J and a are roughly estimated 
at 250 ± 15 K and 0.35 ± 0.05, respectively. 

The competing-J model seems to explain the C m of CuGe0 3 . However, the magnetic 
susceptibility in Ref. ([!]) is not in agreement with the theoretical one with J and a estimated 
in the present study. On the other hand, Riera and Dobry (J = 160 K and a = 0.36)1 and 
Castilla, Chakravarty, and Emery (J = 150 K and a = 0.24)0 have reported different J 
and a, and the theoretical magnetic susceptibility curves with their parameters agree with 
the observed magnetic susceptibility. However, as is shown by the dashed and dotted curves 
in Fig. 5, the corresponding theoretical magnetic specific-heat curves do not reproduce the 
measured one. Here we set C m ^ max by the value of the observed datum at 90 K. Then 
we conclude that we cannot find out the parameters which explain the observed magnetic 
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specific heat and the observed magnetic susceptibility simultaneously. This is a problem in 
the competing- J model. 

Our a (0.35) is close to that of Riera and Dobry (0.36) and seems to be larger than a c . 
This indicates that the energy gap remains finite even above T$p. The energy gap above 
Tsp may be as large as 0.015J which was estimated by Riera and Dobry! when a = 0.36. 
However the finite gap in the uniform state has not been observed by any experiments. This 
is the second problem in this model. 

The competing-J model is an interesting idea and is probably applicable to CuGeO"3. 
However, this model cannot explain quantitatively all the experimental results. The prob- 
lems may be solved by considering (i) interchain exchange interactions, (ii) temperature 
dependence of the exchange interactions, and (iii) structural fluctuations of the lattice dimer- 
ization. 

These suggestions are due to the following reasons; 

(i) In CuGeOs, the strength of the interchain interaction along the b axis is about 10 
% of that along the c axis,i and it cannot be neglected. Riera and Koval suggested that 
the gap did not open in the uniform state even in the case of a = 0.36 when the interchain 
interaction exists, because the interchain interaction changes the critical value a c .i! 

(ii) Because the lattice constant of CuGeOs changes strongly with decreasing 
temperature,l'll'0 J might depend on the temperature. It leads to a substantial effect on 
the temperature dependence of the magnetic susceptibility and the magnetic specific heat. 
When temperature dependence of J is taken into account, both the observed magnetic spe- 
cific heat and the magnetic susceptibility might be explained simultaneously in terms of the 
competing-J model. 

(iii) Pouget et aM observed the structural fluctuations of the lattice dimerization whose 
inverse correlation length is proportional to (T — Tgp) 1 ^ 2 below about 40 K by X-ray diffuse 
scattering. Thus, the fluctuations may affect the magnetic specific heat and the susceptibility 
through a strong spin-lattice coupling. 

Recently, Nakai et aZ.il obtained the magnetic specific heat up to 650 K by means of the 
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birefringence. The differential of birefringence with respect to temperature gives a magnetic 
specific heat. The temperature dependence of the magnetic specific heat shows a broad 
maximum around 80 K. Their result agrees well with ours. 

5. CONCLUSION 

We studied the quasi-elastic Raman scattering in CuGeOs, which is caused by spin- 
fluctuations. The quasi-elastic Raman scattering is observed only in the (c, c) polariza- 
tion. The line shapes of the Stokes and anti-Stokes Raman scattering were well described 
by a Lorentzian curve around cm -1 . This fact indicates that the quasi-elastic scattering 
originates from the fluctuations of energy density in the spin system. The temperature 
dependence of the half width at half maxima, which is proportional to is well de- 
scribed as (T — Tsp) 1 / 2 . This is consistent with the results of the neutron-scattering studies. 
The temperature dependence of k 2 /T 2 , which is proportional to the magnetic specific heat, 
was compared with the theoretical calculation by the competing-J model. The obtained 
J = 250 ± 15 K and a = 0.35 ± 0.05 are different from those obtained in the magnetic sus- 
ceptibility. We cannot find the parameters which explain the specific heat and the magnetic 
susceptibility at the same time. We suggested the reasons of the discrepancies between the 
experimental results and the competing-J model. 
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FIGURES 

Fig. 1. Temperature dependence of the calculated magnetic specific heat for a = 0, 0.10, 0.20, 
0.30, and 0.40 (a), and their normalized specific heat (b). 

Fig. 2. Polarization characteristics of the Raman spectra between -100 and 100 cm -1 in CuGe03 
in the right-angle scattering geometry at 300 K (a). Polarization characteristics of the Raman 
spectra between -100 and 500 cm^ 1 (b). The b(c, c)a and a(c, c)a spectra in (b) are normalized by 
the intensity of the 181-cm _1 A g phonon modes. Ar + denotes a plasma line of the Ar + ion laser 
at -13.4 cm -1 . 

Fig. 3. Typical spectra of the quasi-elastic scattering in CuGeOs at 250 (circles), 120 (squares), 
50 (triangles), and 20 K (diamonds). The solid and dashed curves denote Lorenzian curves and 
Gaussian curves, which are written in Eqs. (6) and (7), respectively. 

Fig. 4. Temperature dependence of T (a) and k 2 (b) above 50 K in CuGeOs. The solid curve in 
(a) is proportional to (T - Tgp) 1 / 2 . 

Fig. 5. Temperature dependence of k 2 /T 2 above 50 K in CuGeOs, and the smoothed data. 
This figure also includes the calculated magnetic-specific-heat curves for a = 0.30 and J = 237 K, 
for J = 261 K and a = 0.40, for the parameters of Riera and Dobry (J=160 K, a = 0.36) (Ref. 
9), and for those of Castilla, Chakravarty, and Emery (J=150 K, a = 0.24) (Ref. 10). All the 
theoretical curves were normalized so that C m ^ max agrees with the observed datum at 90 K (see 
text). 
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